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Properties of an extended periodic Anderson model with the f-c Coulomb interaction [7^ 
is studied as a model for CeCu2Ge2 and related compounds which is cosidered to exhibit a 
sharp valence change under pressure. The problem is treated by extending the slave-boson and 
large- A'^ expansion method to treat the present system. It is shown that, as the f-level et is 
increased relative to the Fermi level, the sharp valence change is caused by the effect of Ufc with 
moderate strength of the order of the bandwith of conduction electrons. The superconducting 
transition temperature Tc due to the valence-fluctuation exchange is estimated on the slave- 
boson fluctuation approximation. In the model with spherical Fermi surface, Tc exhibits sharp 
peak as a function of ef , simulating the effect of pressure, for the d-wave pairing channel. 
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§1. Introduction 

CeCu2Si2 is well known not only as the firstly discovered unconventional superconductor with 
Tc ~ 0.7 K at ambient pressure,^-* but as a remarkable pressure dependence of Tc exhibiting a 
pronounced peak at P ~ 17 GPa.^^ The isostructural compound, CeCu2Ge2 , has a similar phase 
diagram in P-T plane in which one can see the superconductivity appears after the magnetism is 
suppressed by pressure at P 8 GPa.^"^^ These two compounds have similar physical properties, 
if the scale of the pressure for CeCu2Si2 is shifted by 7.6 GPa.^^ The mechanism of such prominent 
enhancement of has not been understood so far. A purpose of this paper is to discuss it from 
the view point that such phenomena is based on the sharp valence change of Ce ion. 

Apart from the pronounced peak of Tc, remarkable properties under pressure of these compounds 
are as foUows:'^^ 1) The residual resistivity po also exhibits a peak at around the pressure where 

* Present address: NEC Corporation, 4-1-1 Miyazaki, Miyamae-ku, Kawasaki, Kanagawa 216-8555, Japan. 
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Tf. takes the maximum. 2) The coefficient A of the T^-term of resistivity rapidly decreases by 
about two orders of magnitude at around the same pressure as Tc and po exhibit peaks. These 
phenomena may be attributed to a rapid change of the valence of Ce ion/^ while the origin of the 
superconductivity around magneic quantum critical point (QCP) is considered to be induced by the 
enhanced spin fluctuations.^^ It is not only becasue the position of the sharp peak of Tc is located 
far away from the pressure corresponding to QCP, but the valence of Ce ion seems to decrease 
rapidly there. The latter fact is supported by the rapid decrease of the coefficient A by about 10~^ 
times, which implies that the Sommerfeld constant 7 decreases by more than 10~^ times suggesting 
that the system is changed from the Kondo regime to the valence fluctuation (VF) regime rapidly. 
This is also supplemented by the rapid decrease of the so-called Kadowaki- Woods ratio by 25 
times there from the value of heavy electrons to that of the conventional transition metals.^' 

Recently, we have proposed that the physics behind these intriguing behaviors should be the 
enhanced valence fluctuations of Ce ion, and given a consistent explanation on the basis of a phe- 
nomenological model for the valence fluctuations.^^ In the present paper we develope its microscopic 
theory on the basis of an extended periodic Anderson model with repulsion between f- and con- 
duction electrons. Without the valence change does not occur so drastically. On the other 
hand, it is shown on the mean-field approximation of slave bosons that the rapid valence change 
can be realized for moderate strength of U{c of the order of the bandwidth of conduction electrons if 
the level of f-electron ef is tuned relative to the Fermi level. It is also shown by taking into account 
the fluctuations of slave bosons beyond the mean-field values that the superconductivity can be 
induced in the d-wave channel associated with the rapid valence change provided that the sherical 
condcution band is adopted. 

In §2, we introduce the extended periodic Anderson model and briefly review the relevant works 
performed so far. In §3, we present a formalism for obtaining mean-field solutions in the slave-boson 
approach and taking into account the fluctuations around it. The results of mean-field solutions 
of slave bosons and superconducting transition temperature are presented in §4. The main resutls 
of this paper are as follows: 1) Sharp valence change is caused by the effect of C/fc with moderate 
strength of the order of the band with of conduction electrons, when the f-level ef is tuned as a 
mimic of the effect of the pressure. 2) The superconducting state is induced by the process of 
exchanging the slave-boson fiuctuations for the values of ef at which the sharp valence change 
occurs. 3) The symmetry of so induced superconducting state is d-wave if the spherical model is 
adopted for conduction electrons. The discussions are given in §5. 

§2. Extended Periodic Anderson Model w^ith f-c Coulomb Interaction 

The effects of pressure appear through changes of the parameters characterizing the physics, 
the Kondo coupling J, the hybridization matrix element V, and the f-level Cf. However, these 
variations of parameters themselves cannot afford causing the drastic change of valence of Ce-ion 
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from that of the Kondo regime to the VF regime, which leads to the divergent increase of the 
valnece susceptibihty, while the valence increases prominently as the strength of V is increased 
or the level ef is increased appraoching the Fermi level. Here we take into account the effect 
of the short range Coulomb repulsion between f electron and conduction electrons, because the 
repulsion may promote the charge transfer between both type of electrnic states in general when 
it is combined with the level shift of Cf. Namely, we extend the conventional periodic Anderson 
model (PAM) as follows: 

kcr ka i 

+yJ2(iafka + h.C.) + C/fc E (2-1) 

ka iaa' 

where the conventional notations for PAM are used except for Ufc, the f-c Coulomb repulsion. 

The effect of U{c on the valence of Ce-ion has been discussed in a number of contexts and 
models. First, with regard to the impurity Anderson model, the effect of U^c was discussed in 
relation to optical experiments, valence-band photoemission (PES) and bremsstrahlung isochromat 
spectroscopy (BIS).^^"^^^ Costi and Hewson studied the effects of U{c by a numerical renormalization 
group (NRG) approach. They concluded that in the Kondo regimes Ufc can be absorbed to 
other parameters {V, ef and U{f), and the same set of renormalized parameters on f/fc = is 
consistent with both the valence-band PES and BIS spectra and the thermodynamic properties. 
The effect of increasing C/pc is to increase the Kondo temperature, thus to decrease the number of 
f-electrons nf per ion mored rapidly. Indeed, a related model has been investigated on the basis of 
NRG approach by Takayama and Sakai, and it turned out that nf rapidly decreases at ef « Ep — Uic, 
E-p being the Fermi level, as Cf is increased so as to approach E-p}'^^ This result implies that the 
rapid valence change occurs where Cf + J7fc, the energy of state, and £'f, the energy of the state 
7*^+ (extra conduction electron on the Fermi level), are nearly degenerate giving rise to enhanced 
valence fluctuations. Then, C/fc has a considerable effect on the valence fluctuation, although it 
does not cause a valence instability in the impurity Anderson model. 

On the other hand, it is not easy to examine the effects of J7fc in the PAM due to the lattice 
effect. In the case of impurity model, the conduction band plays a role of the electron bath so that 
the chemical potential ^ is essentially fixed. However, in the case of lattice model, ^ is considerably 
affected by Ufc itself so that we have to treat the problem in a self-consistent fashion. In particular, 
such a self-consistent treatment is indispensable when the valence state of Cc begins to leave from 
the Kondo regime to VF one under the high pressure. There exist some studies of the extended PAM 
with Uic within Hartrcc-Fock like approximations. Even in this simple level of approximation, 
C/fc is responsible for the rapid change of the number of f-electrons as the level ef is tuned, leading 
to the first order transition in the large f/fc region. Recently, we have investigated the extended 
PAM with f/fc by the variational Monte Carlo method on the extended Gutzwiller variational wave 
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function.^^) It has been shown there that the valence of f-electrons decreases rapidly as the level 
of f-electrons ef is increased, if U^c is moderately large comparable to half of the bandwidth of 
conduction electrons. 

The similar effects of the d-p Coulomb interaction in the so-called d-p model have been inves- 
tigated as a possible charge fluctuation mechanism of the high Tc superconductor.^^' On the 
other hand, the physics discussed in the present paper is rather different from "valence fluctuation 
mechanism" proposed for heavy electron superconductors in Ref. 24. 

On the basis of these considerations, we investigate the effect of U^c in the extended PAM as 
a possible origin of the salient phenomena observed in CeCu2Ge2 under high pressures, partic- 
ularly a possible mechanism of the enhancement of the superconducting transition temperature 
which has not been clarified so far. We study the extended PAM with Ufc by the slave-boson and 
large- A?^ expansion approach,^^"^^) which has been shown to be effective for studying thermody- 
namics properties, superconducting transition temperature, transport properties, etc. of the usual 
PAM.^^"^°) Following these methods, we extend it so as to include the f-c Coulomb interaction 
within a Gaussian fluctuation approximation and study its effects on the superconducting transition 
temperature. 

§3. Formulation 

In this section we summarize a formulation for calculations using the slave-boson large- A ex- 
pansion technique. We generalize the Hamiltonian eq. (2.1) to the case where the A-hold 
degeneracy exists , including both spin and orbital degrees of freedom. Then, we introduce the 
slave-boson operators to exclude the doublely occupied state by f-electrons assuming f/ff = oo. 
Thus, the Hamiltonian is expressed as follows: 

^ = {^kci^Ckm + ^ifijkrr) + {Amfk+qmb]j + h.C.) + Uf^ ^ fufiAmCim, (3.1) 

km ^ kqm il'^ 

where 6's are slave-boson operators describing f ° state and the constraint 

Qi = T.fimfim + b\bi = l (3.2) 
m 

is required at each site in order to maintain an equivalence of the truncated Hilbert space and the 
original one. To generate a 1/A-expansion we rewrite variables as follows: 



Qi 


^ qoN 


b 


bVN 


V 


v/Vn 




^ f/fc/\/A 



(3.3) 



Hereafter, we use the radial gauge following Ref. 25. (Although the radial and the Cartesian 
gauge formulation are ultimately equivalent, spurious infrared divergences do not appear in the 
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radial gauge approach.^''')) We perform a local gauge transformation, 6i(r) = /9j(r)e'^'('^), fimij) = 
//^(r)e'^*('^) and A^(r) = Aj + ^i(r). Then, the partition function is given by 



where 



Z = j Vicc^fppX)eM-S), 



(3.4) 



S= dr 
Jo 



(3.5) 



In the above expressions, we have rewritten variables /' and A' as / and A, respectively, and 
neglected the Jacobian factor Hir following Ref. 25. 

By introducing two kinds of Stratonovich-Hubbard fields (p^ and (p'^ for Ufc, we can perform the 
functional integral over the Fermion fields in eq. (3.4) and eq. (3.5). Then, the partition function 
eq. (3.4) can be transformed as follows (See Appendix A A for a detail of derivation): 



Z = y"p(pA¥?V')exp(-5), 



(3.6) 



N 



S = - NTTlnA + i^=T^ V p(-k)A(k' - k)p(k) 



kk' 



iqoNVK J drA(0;r) - i^T^ (k)<^^(-k), 



(3.7) 



where we have introduced abbreviation k = [k, icOn), etc. and matrix A is defined as 



^kk' 



i-i^n + efe)<5uk' + ^^^'(k - k') 



p(k-k') 



T^P*(k - k') {-iun + ef)4k' + ^iA(k - k') + |^<^^(k - k') 



At the level of mean-field approximation, the action So can be written as 

So = -A^IVlnio + - Qo^/N[yX - ^^^^(k)^^(-k) 



(ioJn — ei 



(p I P 

iVL 2VnE Ni^ 



(3.8) 



(3.9) 



(3.10) 
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where we have approximated A, p and by their uniform mean field value as 

A(q) = ^A5q 



(3.11) 



and 



^Okk' = 



I 1 t'tc -f 



V 



- iUn + €i + 



+ 



3kk' 



= -Gn 



(3.12) 



Here Go is nothing but a matrix of Green's functions for the renormalized band. In this mean- 
field approximation, quasiparticles acquire the heavy effective mass through the strong f-f Coulomb 
interactions, Uff = oo, and deep f- level Cf far below the Fermi level. Prom the optimum conditions 
for mean fields, dSo/dp = 0, dSo/dX = and dSo/dip = 0, we obtain the following set of self- 
consistent equations: 

iA T ^ y2 



Qo - ( 



= --E 



iVL^(ia;„-e-fc)(ia;-e-f)-y2' 



P 



T 



lUJr, 



ek 



iVL V(i^n-efe)(ia;-6-f)-F2' 



efc 



<p 



where we have introduced abbreviations as 



-y- 

Nj^^ {iu;n-ek){iio-ei)-V'^' 
-J: 



"wn - efe)(it^ - Cf) - ' 



ek = ek + Uicni, 

€f = ef + iA/ a/ZVl + [/fcflc, 



(3.13a) 
(3.13b) 
(3.13c) 
(3.13d) 

(3.14a) 
(3.14b) 

(3.14c) 



where nf = {nf)MFi the f-electron numbers per site and "spin", and Uc = (?^c)MF^ the number of 
conduction electrons per site and "spin". It is noted that both the left-hand sides of (3.13b) and 
(3.13c) are equal to rlf and that of (3.13d) is equal to fie- Thus, the following relations hold: 



qo 



2VNl 
fie- 



nf 



In terms of these quantities, the matrix Green function Gq, eq. (3.12), can be expressed as 



Go 



1 



{iujn - efc)(iw - ef) - V'^ 



luj — ef 
V ioj 



V 



(3.15) 
(3.16) 

(3.17) 



6 



'-^0 



(3.18) 



Next, we present the formalism treating the Gaussian fluctuations around their men-field values. 
For this end, we express the boson fields as sums of the mean-field value (indicated by "bar") and 
the fluctuations around it (indicated by "tilde") as 

A(q) = iA<5q + A(q) 

p(q) = TP'^q + p(q) • 

V?(q) = ^(p5^ + (^(q) 

Then, we expand S, (3.7), with respect to the fluctuations up to the second order as 



(3.19) 



S = Sq + S, 

where Sq is given by (3.10) and the fluctuation part S is given as follows: 

~S = -NTrlnA+^TY,{Xp{-k)p{k) + 2p~X{-k)p{k)} - ^T^0'{k)0%-k) 



(3.20) 



(3.21) 



where A= 1 — GqM with M = A — Aq, the explicit form of which is 

"^(p\k-\^) T-J - 

T V 



^p*(k-k') r^A(k-k') + r^<^^(k-k') 



r^p(k-k') 



The first term of eq. (3.21) can be expanded as 

oo 

- Tr In i = - Tr ln(l - GqM) = ^ Tr(GoM)". 



(3.22) 



(3.23) 



n=l 



We truncate this expansion at n = 2 obtaining a Gaussian form. Then, we can write the partition 
function of the Gaussian fluctuation part as (See Appendix B) 

Z = j V{p\(p^(p'')ex.Y>[-S], (3.24) 

S = iVT ^(p(-k), iA(-k), ^^(-k), (-k)) 5k (p(k), iA(k), r (k), ^\k))\ (3.25) 

k 

where the superscript t means the transpose is taken, and the symmetric matrix Sk is given as 

Spp{k) S'pA(k) Sp^.{k) Sp^t{k) 

S\p{k) S'AA(k) 5'A<^<=(k) 5'A^f(k) 

S^-p{k) S'^cA(k) S'^c^c(k) 5^c^f(k) 

_S^tp{k) S^fy^{k) S'^f^c(k) 5^f^f(k) 



5k 



(3.26) 



iA 

/JVl 



+ yW(k) + y2n3(k) 

-k 



+ VYi2{k) 

^nf(k) 



%yn2(k) 

C/fcTTf 



n*(k) 



%yni(k) 
s-n^f(k) 



^n^k) -^ + --^w\k) 
fn^(k) 



(3.27) 
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where the definition of polarization functions II's are given in Appendix B, and the symmetric off- 
giagonal components have abbriviated. The inverse matrix of gives the matrix Green Functions 
for bose fields 

D^fs{k;T) = -(T^a(fe,r)/?t(fc,0)), (3.28) 

where a, P are p, A, ip^, ip^, as 

£)(k) = -S-^k). (3.29) 

In Fig. 1, we show the various possible interaction terms included in eq. (3.21) by Feynman dia- 
grams. 



Fig. 1. Feynman diagrams of interaction terms of eq. (3.21). (a) represents the hybridization process, (c), (d) and 
(g) are new terms which appear through the f-c Coulomb interaction (e) and (f) represent the second term of 
eq. (3.21), and (g) corresponds to the last term. 



§4. Results 

4-1 Mean-Field Solutions 

First, we present the solutions within the mean-field approximation. Self-consistent equations 
(3.13) can be rearranged as follows:^^) 

q^-^_^f = A^y[l±^^JlZlh^^]f^E^) (4.1b) 

j^E{/(^fc)+^(^fc)> = ^^+^^ (4.1c) 

k 

where f{x) is the Fermi distribution function and are the quasi particle energies 

E% = hei + e-fc ± ^(6-f-6-fe)2 + 4F2]. (4.2) 
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A set of equations (4.1) is equivalent to that for the case of conventional PAM without Ufc except 
for the point that and Cf include the term arising from Ufc as seen in (3.14). We solve them for a 
three-dimensional model with the dispersion for the conduction-electrons approximated by that 
of free electron, e^. = fc^/2m - D, where the bottom of the conduction band is set as —D. Hereafter 
we use D as the unit of energy. The bare mass, m, of conduction electrons is chosen such that the 
integration of, po(e), the density of state per "spin" of conduction electrons with respect to e from 
—D to D is equal to 1: 



In Fig. 2, we show the results of rif as a function of ef for several values of Ufc/D. Here we 
set the hybridization as V =0.5D and total electron number per "spin" as n = rlf -|- nc = 0.875. 



physical picture that the rapid valence change occurs at ef ~ Ep — Ufc where the energies of Z'^- 
and /^-state are degenerate leading to enhanced valence fluctuations. For much larger values of Ufc 
or smaller values of V than those presented in Fig. 2, there occurs a first-order like discontinuous 
transitions although they are not shown in Fig. 2. In the mean- field level of approximation for 
slave boson Hamiltonian, our treatment of the effect of C/fc is just like that in the Hartrcc-Fock 
approximation. It is noted that the valence change occurs more sharply if we estimate it in 
much more proper approximation on the extended Gutzwillcr variational wave function. In this 
sense, the sharpness of the valence change may be underestimated by the present treatment. 



Fig. 2. m, f-electron number per site and "spin", as a function of ef, f-electron level, in three-dimensinal model 
where the free electron dispersion = /2m — D is adopted for the conduction electrons. 

Using the mean-field solution we calculate the density of states p(//) of quasiparticles at the Fermi 




(4.3) 



These results are consistent with those of previous works. "^^^ These results are consistent with a 
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energy. It is related to the Sommerfeld constant 7 of the specific heat as 

7 = y i^p()U)iV, (4.4) 

where N is the degeneracy of the quasiparticles. The so-called Kondo temperature, or the charac- 
teristic temperature, of the present model, is defined simply as (ef — n). In Fig. 3, we can see that 
the relation Tk oc 7"^ holds. It is remarked that the relation Tk vs. 7 with different value of Ufc 
are lying on a line exhibiting a kind of universality. It is suggested that the effect of C/fc can be 
absorbed into the other parameters of the conventional PAM at least in the mean-field level, like 
in the case of the single impurity Anderson model. However, in the lattice case, the first-order like 
tansition occurs through the effects of Ufc even in the mean-field approximation. This is considered 
to be a distinctive effect of U^c in the PAM. 



Fig. 3. Sommerfeld constant 7 vs Kondo temperature Tk, defined by Tk = et — //. The parameters used are the 
same as those in Fig. 2. Results for different values of Uic/D lie on the same line, exhibiting a universal behavior. 



4-2 Superconducting Transition Temperature 

Next, we discuss a possible type of the superconducting gap near the region where the rapid 
valence change occurs owing to the f-c Coulomb interaction in the extended periodic Anderson 
model. In the conventional PAM, problems of determining the superconducting transition temper- 
ature have been studied by several authors within the slave-boson and 1/A^-expansion method. ^^'^^^ 
Following their method, we calculate the superconducting transition temperature in the weak cou- 
pling theory in the present framework. The transition temperature is obtained from the integral 
equation of the particle-particle scattering amplitude, F, for two quasiparticles with opposite mo- 
mentum near the Fermi surface. (See Fig. 4.) The quasiparticle operators can be represented by a 
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Fig. 4. Fcynman diagram representing of the integral equations for the scattering ampUtude of two quasiparticles 
rfrom (fc,-fe) to (fc',-fe'). 



Fig. 5. Feynman diagram for the irreducible vertex part P'"' to leading order in 1/A'^. The explicit form of the 
fluctuation propagators D is given by eq. (3.29) and its interaction vertices g with quasiparticle is shown in Fig. 1 



unitary transformation in terms of f- and conduction electron operators as 













ffe 















(4.5) 



where 7*^^'' corresponds to the eigen values , eq. (4.2), respectively. The scattering amplitude T 
is obtained from two-quasiparticle correlation function, 

(T.7L"' M-^Wt' <(^3)7L-'^(r4)), 

K ,m —ft ,m —K,m K,m 

by removing the external legs. 

To leading order in 1/N, irreducible vertex part T^^^ includes only a single-boson exchange process 
as shown in Fig. 5. With the use of the relations, eq. (4.5), eq. (3.29), and interaction vertices shown 
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in Fig. 1, analytic expression of r(°) is given as follows: 



, 2 2 
+ U V 



pp 



(4.6) 



In the weak coupling limit, in which the external momenta are set on the Fermi surface, i.e. 

— kp, and the static limit, a; — ^ 0, is taken in the boson propagator, the linearized gap 
equation represented by Fig. 4 can be solved by a conventional method. Namely, the scattering 



Fig. 6. To for d-wave channel and rif, f-electron number per site and "spin", as a function of ef. 



amplitude is decomposed into the Legendre polynomial as 

oo 



(4.7) 



The scattering amplitude T( corresponding to the channel with relative angular momentum £ is 
given by 



re 



l+p(/x)rf ln(^)' 



where F^'^^'s are related to r('^)(fe • fe )'s by the formula 



Vf^ = \ d(cos0)r(°)(e)Pi(cos^). 
2 J -I 



Then the transition temperature for ^-wave channel is given by 

rc = rKexp ^ 



Lp(M)r 



(0) 



(4.8) 



(4.9) 



(4.10) 



Here it is noted that the energy cut-off, corresponding to the Debye frequency, is given by Tk = 
ef — ji, the bandwidth of quasiparticles. 
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The transition temperature so calculated can exist only for the d-wave {£ = 2) channel as 
far as the channels, £ = 0, 1 and 2, are concerned. In Fig. 6 we show as a function of ef , the 
f-electron level relative to the Fermi level, for seeral values of C/fc, the repulsion between f- and 
conduction electrons. Parameters adopted are the same as those in Fig. 2. There exists a sharp 
peak of Tc at around ef where nf starts to show a rapid decrease. Its tendency becomes more drastic 
as Uic increases making the valence change sharper. In the region where the f-electron number is 
decreased enough, T^. is strongly suppressed. We have investigated which term in eq. (4.6) plays 
important role for pairing interaction. It turned out that the major part of T^^^ is induced by the 
scattering process (f,f)— >-(f,c) or (f,c)— *-(f,f), in which the valence of f-electrons is changed directly. 
(See fig. 7.) 



Fig. 7. r"^"^ (g) as a function of momentum transfer q for each scattering channel included in eq. (4.6). For example, 
"fffF' means the term of eq. (4.6). The others are represented similarly. 



In the Kondo limit where nf 1/2, the spin fluctuations are believed to play the most important 
role for the occurrence of superconductivity. In such a region we have to take into account the 
higher order term beyond 1/N to discuss the instability to the superconducting state, since the 
spin-fluctuation contribution to the effective interaction appears only beyond at the order (1/A^)^.^°) 
However, the present approach of the order of l/N is still expected to work in the region where 
the valence fluctuations play an important role. 

§5. Conclusions and Discussions 

We have developed a microscopic theory to support the idea that the rapid change of f-valence 
of Ce ion is the origin of anomalous behaviors of CeCu2Ge2 under the pressure P ~ 17 GPa, such 
as the coincidence of peaks of Tc and po, and the rapid decrease of the coefficient A of T^-term of 
the resistivity. To this end, we analyzed the properties of an extended periodic Anderson model 
with the f-c Coulomb repulsion C/fc by the slave-boson and 1/A^-expansion method. It turned out 
that such a model indeed contains the ingredient making rapid change of nf , f-electron number per 
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site, possible. It was also shown that the superconducting transition temperature of d-wave pairing 
is sharply enhanced just before the rapid decrease of nf as ef is increased. It turned out that the 
scattering process of (f,f)— >-(f,c) or (f,c)— >-(f,f), in which the f-electron number is changed directly, 
plays a major part in Cooper pair formation. Cooper pairing of channels other than the d-wave 
cannot be possible so long as the pairing interaction discussed in the present paper is concerned. 

A few remarks should be made concerning the following points. First, although we did not men- 
tion the relation between the coefficient A of T^-term of the resistivity and the Sommerfeld constant 
7 within the present approach, we have verified that the Kadowaki- Woods relation is reproduced 
in the heavy electron limit.^^^ However, the crossover of the universal ratio A/j"^ between classes 
of the heavy electron and of the transition metal was not reproduced. We have cognizance of the 
necessity of calculations in which the dynamical effect of quasiparticle selfenergy is fully taken into 
account on the order of 1/N. However, we have to left such calculations for a future study because 
we need much more considerable amount of study for performing its program. 

Second, there exist several points of view different from ours for explaining the anomalous proper- 
ties of CeCu2Ge2 . One of them is to attribute its origin to the orbital fluctuations in the multiband 
periodic Anderson model, in which the broad bandwidth under high pressure is expected to change 
the degeneracy of the f-electron state. Indeed this has been proposed in ref. 4 paying attention 
to the fact that, at the pressure corresponding to the peak of Tc and po, the two temperatures 
T™^ and T2^^ at which the resistivity takes peak merges with each other, indicating the Kondo 
temperature becomes of the order of the crystal field (CF) splitting. The orbital fluctuation mech- 
anism has also been proposed as a possible mechanism for explaining the phenomena observed 
in CeCu2Ge2, and discussing that the orbital fluctuation is enhanced where Tk is comparable to 
CF splitting.^^) However, it is not clear whether all the anomalies observed in CeCu2Ge2 can be 
explained by those mechanisms. 

Third, the origin of superconductivity around the QCP is believed to be due to the enhanced 
spin fluctuations, while we have neglected here its effect because we are interested in the physics in 
the region far from the QCP. Therefore, the primary origin of the superconducting transition may 
change to the spin-fluctuation mechanism as QCP is approached. In order to discuss this crossover 
in the present model, we have to study the problem up to (l/AT)^ which is beyond scope of this 
paper. 

Last, we have discussed the superconducting transition temperature in the weak coupling theory. 
However, it is suggested that the frequency dependence of the selfenergy is important as can be 
seen in the discussions concerning the Kadowaki- Woods relation.^^^ So, the study in the strong 
coupling theory of superconductivity is desired for more solid conclusion. Nevertheless, the result 
does not seem to be modified qualitatively, judging from our past experience in the problem of 
Cooper pair formation in isotropic two-dimensional Fermions with short range repulsion.^^^ 
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Appendix A 

Here we discuss how the two kinds of boson fields are introduced in relation to the f-c Coulomb 
interaction. The f-c Coulomb interaction term in eq. (3.1) can be rewritten as 



(Al) 



With the use of the identity 



exp 



AN 



f dr{(-ir/i(r)+E(4a + <a))' 

-(-a(r)+EK'a-<a))'} 



= Const. 



(A-2) 



and eq. (A-2), the term including C/fc in eq. (3.4) can be expressed as follows: 



V{ri C) exp 
= j P(cct//tryC)exp 



a/3 

'dr{^(r?f(r) + C^(r)) 



+ ^(Mr) - a(r)) E^L + (i^^W + C.(r)) E<a)} 

a 

= 1^(00^/^^=) exp -l^dr{ 



(A-3) 
(A-4) 



+ ^(¥'?(r)E"L + ^'(^)E<a))} 



(A-5) 



In the last line, bose fields, 77 and have been transformed to the new ones, (fl = N{irii + 
and ipj^ = N{irii — C,i)- After carrying out the functional integral over the Fermion fields, we obtain 
eqs. (3.6) and (3.7). 
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Appendix B 

Here we derive formal expressions of the Gussian fluctuation part of the action, eq. (3.21), and 
symmetry property of the polarization functions, eq. (3.27). Of the terms in the part of action, 
eq. (3.21), 

-Trlni = -TVln(l-GoM) = ^Tr(GoA^)", 

n 

those with n = 2 give contributions of the Gaussian fluctuation. Namely, the relevant term is 
expressed as 

V(GoM)2 = ^tr{Go(k)Mkk+qGo(k + q)Mk+qk}- (B-1) 

^ kq 

The explicit form of GqM is given as 

GoM(k,k + q)= (B-2) 

T r %G8^(k)^f (-q) + VGl'{k)p{-q) VG^{k)p{q) + Grf(k)(iA(-q) + %<^^(-q)) ' 
VNZ [ S.Gfc(k)<^f (-q) + FG«(k)p(-q) VG'^{k)p{-q) + G«(k)(iA(-q) + ^^'^(-q)) 

Then, the diagonal components of the product of four matrices in eq. (BT) are written as follows: 



[G'oM(k, k + q)G'oM(k + q, k)](i i) = (B-3) 
^ f{^G^^(k)^f (-q) + FG^f (k)p(-q)}{^G^^(k + q)^f (q) + FG^^(k + q)p(q)} 



+{yGS^(k)p(-q) + GS^(k)(iA(-q) + ^^'^{-q))}{^Gi^{k + q)/(q) + FG«(k + q)p(q)}j ^ 



and 

[GoAf(k,k + q)GoM(k + q,k)](2 2) = (B-4) 
^ (^{^Gfc(k)^^(-q) + FG«(k)p(-q)} 

{VGl%k + q)p(q) + Gfik + q)(iA(q) + ^^^(q))} 
+{VGl%k)p{-cD + G«(k)(iA(-q) + ^^^(-q))} 
{yG^o'^(k + q)p(q) + Gl'{k + q)(iA(q) + ^^^(q))} j • 

Substituting these expressions into eq. (B-1) and making rather tedious rearrangements, we obtain 
eq. (3.25) with the matrix S whose components are deflned by eqs. (3.26) and (3.27). Polarization 
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function II's in eq. (3.27) are defined as follows: 



(A;,icj) 


= n(*^'ff)(fc,icj) 






n^^(A:,iu;) 


= n(^f'^f)(fc,ia;) 






Ili{k, iuj) 


= i{n(^^'^f)(fc,ic^; 




iuj)} 




= i{n(ff'^f)(fc,iw) 




ico)} 


U3{k, iuj) 






ia;)} 


as 









(B-5) 



(B-6) 



where a,(3,j and (5 stand for c or f. It is noted that n^"^'''''') satisfies the following symmetry 
relations: 

n(^^'^f)(A;,ia;) = n(^f'^'=)(fc, -io;) 

n(ff'^f)(fe,ia;) = n(^f'ff)(A;,-ia;) (B-7) 

n('=^''f)(A;,ia;) = n(ff'^^)(fe, -iw). 
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